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BI-LEVEL INTEGRATED SYSTEM SYNTHESIS (BLISS) EOR CONCURRENT 
AND DISTRIBUTED PROCESSING 

Jaroslaw Sobieszczanski-Sobieski + . Troy D. Alius*. Matthew Phillipsi, and Robert Sandusky* 


Abstract 

The paper introduces a new version of the Bi- Level 
Integrated System Synthesis (BLISS) methods intended 
for optimization of engineering systems conducted by 
distributed specialty groups working concurrently and 
using a multiprocessor computing environment. The 
method decomposes the overall optimization task into 
subtasks associated with disciplines or subsystems 
where the local design variables are numerous and a 
single, system-level optimization whose design vari- 
ables are relatively few. The subtasks are fully autono- 
mous as to their inner operations and decision making. 
Their purpose is to eliminate the local design variables 
and generate a w ide spectrum of feasible designs whose 
behavior is represented by Response Surfaces to be ac- 
cessed by a system-level optimization. It is shown that, 
if the problem is convex, the solution of the decom- 
posed problem is the same as that obtained without 
decomposition. A simplified example of an aircraft 
design shows the method working as intended. I he 
paper includes a discussion of the method merits and 
demerits and recommendations for further research. 

Nomenclature 

AR aspect ratio 

B B “black box", another term for a module 
c coupling equality constraints 

g behavior constraints local to a BB 

h equality constraints tantamount to solution of 

analysis 

ESF engine scale factor 

h cruise altitude 

L/D lift to drag ratio 

Q system level design variables, j Xsh | Y* | w[, 

a subset of Z 
R flight range 

RS Response Surface, SRS - sheaf of RS 

SCF specific fuel consumption 

Sref wing reference area. 

t/c airfoil depth 

ts the cross-sectional dimensions 


f NASA Langley Research Center. MS 240, Hampton, VA 23681. 
♦The George Washington University, the LaRC campus. 


T throttle setting 

U local design variables, JXIoc | Y A J. a subset 

of Z 

U & L upper and lower bounds 

Lht horizontal tail location coordinate 

M Mach number 

w weighting factor in suboptimization objective 

function 
W weight 

Wt total weight, also denoted by TOGW tor take- 
off gross weight 

Xloc design variables local to a BB 
Xsh shared design variables affecting directly two 

or more Black Boxes (modules), BB 
y* behavior variables input into a BB from other 
BB s 

y A behavior variables output from a BB, some 

elements of Y A are designated to be Y* 

Y A S a particular data item selected in a particular 

BB output to be the system objective 
Z - | Xsh | Xloc | Y* | Y A } as a vector of variables' in 
a not decomposed, combined analysis and 
optimization problem 
X taper ratio 

A sweep angle 

0 effective wing area change due to twist 

Superscripts: 
a approximate 

A output from a BB 

* input to a BB generated by the system 

optimizer 

Subscripts: 

0 optimal 

E engine 

F fuel 

HT horizontal tail 

s system objective 

T Take-off-Gross-Weight, TOGW 

w ? wing 


1 Uxamples in aircraft design: Xsh wing aspect ratio and sweep 
angle: Xloc wing cover panel thickness and a composite plv 
oriental ion angle: Y* and Y A an elastic deformation that alters the 
wing aerodynamic shape, and aerodynamic loads that cause the 
wing deformation. 
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0. Introduction 

The relentless drive of computer technology toward 
ever higher computing speed (Floating Point Opera- 
tions Per Second, FLOPS) has enabled solution of 
large computational problems encountered in design of 
complex engineering systems, e.g., an aircraft, in a 
small fraction of time that was required just a lew years 
ago. And yet. it is still not practical to solve such prob- 
lems in one monolithic calculation because the number 
of design variables is large, the analysis is high fidelity 
and it consists of a system of coupled codes. Such 
monolithic calculation, perhaps, may never become a 
routine practice regardless of the computing speed 
available because design of complex engineering sys- 
tems requires that all sorts of information ranging from 
experimental to computational be synthesized by hu- 
man judgment. Consequently, it is a common practice 
to conduct a design process by collaboration of 
autonomous groups of specialists that retain control 
over their domains of expertise, and work concurrently 
to compress the project elapsed time. 

Motivated by this, several methods have emerged 
to enable decomposition of the system optimization 
into a set of smaller tasks aligned with disciplines or 
physical subsystems. Recently, it became evident that 
the above methods are natural for computing w ith large 
numbers of concurrently operating processors (Mas- 
sively Concurrent and Distributed Processing, MCDP), 
a new technology now developing a new computational 
infrastructure having a potential to support engineering 
design at an unprecedented level of effectiveness. 

Optimization by decomposition, initiated by 
Ref 1, became a key component of the Multidiscipli- 
nary Design Optimization- -T 4 , and resulted in a body 
of literature much too large to be cited here in full. A 
few representative references are Refs. 5-15, and 16,17 
(the last two are surveys). Their common thread con- 
sists of suboptimizations at a lower level, all coordi- 
nated by solution of a higher, system-level optimiza- 
tion. Their diversity lies in the means that link the two 
levels, in the use of various approximations, and in the 
definition of the objective function(s) at both levels. As 
to the linking, the typical approaches are the sensitivity 
of optimum to parameters, the system behavior sensi- 
tivity. control over the contributions a suboptimization 
makes to the satisfaction of constraints in another 
suboptimization, and imposition of targets on the 
suboptimization results. 

The approximations commonly used in construct- 
ing surrogate linear or nonlinear models are often based 
on the disciplinary and system behavior derivatives. 


Recently, domain-spanning approximations (DSA) such 
as the Response Surfaces (RS) or Neural Nets (NN) 
constructed with the aid of Design of [Experiments 
(DOE) techniques became popular because they lead to 
repetitive but independent calculations that can be per- 
formed simultaneously to exploit the MCDP technol- 
ogy (e.g.. Refs. 8, 1 1. and 18). 

I he choice of the objective function(s) is funda- 
mentally important. It is also difficult because in a 
suboptimization the objective ought to reflect local 
information combined with the influence of that infor- 
mation on the system performance, and vice versa at the 
system level, while preserving the suboptimization 
autonomy. 

A structural design of the w ing is an example that 
illustrates considerations that underlie the choice of the 
objective in a subsystem optimization. Its output in- 
cludes structural weight and elastic deformation. Con- 
ventionally, the weight would be the objective. How- 
ever, the ultimate objective is a measure of the aircraft 
system performance. That performance, in general, 
benefits from reduction of weight and drag. Because 
increased structural stiffness lowers the elastic wing 
drag, it follows that, for the system benefit, the wing 
structure ought to be designed for an optimal combina- 
tion of low weight and high stiffness. Furthermore, the 
wing structural weight may be lightened by tailoring 
the wing deformation to redistribution of the aerody- 
namic load so as to reduce the root bending moment. 
The resulting load distribution is. in general, different 
from the one that minimizes the drag. 

The original version of BLISS 19,20 addressed the 
above considerations by computing the system sensitiv- 
ity derivatives of the system objective function with 
respect to the subsystem design variables and. then, 
define a subsystem objective function as a sum of 
the subsystem design variables weighted by these 
derivatives. The subsequent versions included the 
use of RS IS and a variant specialized for structural 
optimization- 1 . 

The key new concept in BLISS 2000 reported 
herein is a new formulation for the objective function in 
the subsystem optimization and the use of that optimi- 
zation to control the subsystem output for the benefit of 
the system performance. In the new formulation, the 
subsystem optimization objective is a sum of the sub- 
system outputs, each output weighted by a coefficient 
that is treated as a design variable in the system-level 
optimization. This formulation eliminates the costly 
system sensitivity and optimum sensitivity analyses. It 
also enables representation of the subsystem optima at 
the system-level by the Response Surfaces prepared 
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“off-line” by autonomous, concurrent, and potentially 
distributed operations. It is expected that the specialty 
groups in charge ol these operations work in a MDCP 
environment, being tree to draw on any sources ot 
information and to apply any tools ot their own choice. 

The remainder of the paper defines the BLISS 2000 
algorithm first. It follows with a discussion and inter- 
pretation of the algorithm in qualitative terms, and 
offers a theoretical justification ot the approach by 
demonstrating that the solution obtained by decomposi- 
tion also solves the same problem formulated without 
decomposition. Remarks about the computational cost 
conclude the method description. 

A simplified, conceptual-level design problem ot a 
supersonic business jet illustrates the method and con- 
tributes an empirical basis to a concluding assessment 
of the method limitations, advantages, and recom- 
mended future developments. 

L The BLISS 2000 Algorithm 

The object of optimization in BLISS is a generic 
engineering system whose representative example is an 
aircraft (Fig. 1). The diagram in Fig. 2 shows a few 
major disciplines that form a coupled system in the 
aircraft mathematical model. The model is comprised ol 
modules, also called Black Boxes (BB), coupled by 
data exchanges. 

LI. The Original Problem 

For any optimization-bv-decomposition scheme to 
be valid, the solution it generates should also so lve the 
problem formulated without decomposition . Therefore, 
introduction of the BLISS algorithm begins with tin 
all-in-one formulation that combines analysis and opti- 
mization of a generic system ot which Fig. 2 is an ex- 
ample. It reads 


FIND: Z 0 

MINIMIZE: F(Z) - Y A (Z) a) 

SATISFY: g(Z) <= 0, tor each BB; b) 

h(Z) = 0, for each BB; c) 

c(Z) = Y*-Y A = 0 d) 

ZL <= Z <= ZU; side constraints. e) 
OUTPUT: Z 0 , F 0 A 


In the above, the inequalities g represent the behavior 
constraints local to a BB : , and the equalities h corre- 
spond to the solution of the BB governing equations 
(the BB inner analysis)' . The output-input equalities c 


: Example: structural strength; limit on the aerodynamic pressure 
gradient along the wing airfoil chord. 

'Example: equations of the finite Element Analysis. 


describe the system inner couplings 4 , and in Eq. 1 they 
are present for two reasons: 

- The input Y* to a particular BB is supplied by 
the optimizer instead of being received directly as Y A 
from another BB, where it would originate in a conven- 
tional analysis of a coupled system. Therefore, the 
equalities c are needed to ensure that upon convergence 
the optimizer-generated Y* is the same as the corre- 
sponding Y A . 

- The system analysis needs to be solved. Equating 
the corresponding pairs Y* and Y A is tantamount to 
obtaining a solution. 

The formalism of Eq. 1 is known as the Simulta- 
neous Analysis and Design (SAND). It brings the h and 
c equalities under a single optimization algorithm in 
contrast to the more common optimization techniques 
that eliminate Y A and Y* from Z by solving the h an c 
equalities for each set of the trial values of Xsh and 
Xloc sent out by the optimizer. 

12. BLISS 2000 Optimization by Decomposition 

The BLISS 2000 algorithm separates the BB 
suboptimizations from the system optimization. To 
introduce the algorithm, this section provides a formal 
statement for the optimization at the BB level, fol- 
lowed by definition of the Response Surfaces that link 
the above optimization to the system-level optimiza- 
tion, and concludes with a formal statement lor the 
latter. The section that follows elaborates on a rationale 
for the two-level procedure, its theoretical underpin- 
nings, and the computational cost considerations. 

12 1. BLISS Subsystem (Black Box) Optimiza- 
tion Problem Repeatedly Executed tor a Number o l 
Dispersed Design Points . The design variables Z arc 
divided into the system level variables {Xsh I Y h ) and 
the subsystem level variables {Xloc I Y A ) and the fol- 
lowing optimization problem is solved for each BB. 


GIVEN: Q - {Xsh | Y* | w j 2) 

FIND: U - {Xloc j Y A j a) 

MINIMIZE: «U) = Sw i Y A I *>) 

SATISFY: g(U) <= 0, for each BB; c) 

h(U) = 0, for each BB; d) 

UL <= U <= UU e) 

OUTPUT: Y A o, and Xloc n 0 


The weighting coefficients w, in t, one coefficient 
per element of the output vector Y A , are appended to 
the system-level variables in order to link the BB 


4 Example: the aerodynamic loads that deform a wing should be 
computed for the wing deformed due to those loads. 
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suboptimization to the system optimization. For rea- 
sons to be explained later, the constrained minimum of 
t is not passed to the system optimization; only the 
Y A n values are transferred. 

I he Eq. 2 optimization is executed at a number of 
points dispersed in the Q-space that belongs to the BB. 
The formal Design of Experiments (DOE) techniques 
may aid in forming the dispersal pattern to achieve a 
reasonable coverage of the domain. From the system 
perspective the solution method choice is immaterial 
and that choice need not be the same for all the BBs. 
The specialists in charge of a BB are free to use any 
suitable method, even experiments, or a guess. 

_LL Approximate Model for Optimized Subsystems 

Next, a Response Surface (RS) is fitted (or a Neu- 
ral Net may be trained) to each of the elements of Y A 1S . 
Regarding each RS as a leaf in a sheaf, the resulting 
data base is a Sheaf of RS (SRS). It constitutes an ap- 
proximate model in the Q space of the BB optimized 
using the U coordinates, so that 

Y v, u - Y Ai, 0 (SRS(Q)); and U a « = U a „(SRS(Q)) 
for QL Q <=■ QU (3) 

In the above, the superscript kk a" denotes approxi- 
mate values, and the notation such as Y A<, „ = Y a!, (> 
(SRS(Q)) means that the approximate values are re- 
trieved from precomputed SRS. The bounds on Q, QL 
and QU, are the best estimates, accounting for any side 
constraints and, also, incorporating the move limits 
necessary for the ensuing inter-level iteration. The SRS 
may be called a domain approximation because it cov- 
ers the entire Q space within the above bounds. 

In principle, a similar RS approximation could be 
constructed for each element of Xlocu. However, the 
resulting volume of data to be stored may be so large as 
to make this impractical. If so, Xlocu can be regenerated 
in a manner discussed in Sec. 1.5. 

The computational cost of an RS is proportional to 
the number of points to be evaluated which grows with 
the number of variables. It is desirable, therefore, to 
reduce the number of variables in Q by a condensation 
technique. Such condensation is imperative for those 
variables in Q that represent field quantities, e.g., the 
field of the pressure loads distributed over the wing, 
and the corresponding field of displacements. A con- 
densation of a field data may be accomplished by defin- 
ing the field variable P as 

P = P(p) (4) 


where the number of the parameters p is made as small 
as possible . The p parameters enter the procedures as if 
they were Y* and Y A , each p being represented by its 
own RS. 

When all the BBs are optimized, the functions in 
Eq. 3 for each BB are available to the system optimiza- 
tion that executes next. 

1.4. BLISS System Optimization 

The system-level optimization problem is; 


GIVEN: a set of SRS, one for each BB 5) 

FIND: !Qj - |Xsh | Y* | wj a) 

MINIMIZE: F(Q) - Y A V b) 

SATISFY: c = (Y*,Y A (Q)) = 0; c) 

QL <= Q <= QU d) 

OUTPUT: Qo, F„ e) 


Eq. 5 may be solved by any optimization tech- 
nique, the data for F(Q) - Y Aa « K and c(Q) for any Q 
being retrieved, as per Eq. 3, from the SRS data base. 
The efficiency of the technique is not critical because 
that data retrieval is nearly instantaneous. As in Eq. 1. 
c - 0 in the above is tantamount to a system analysis, 
thus the BLISS algorithm implements, in effect, the 
Simultaneous Analysis and Design method at the sys- 
tem level, but not necessarily at the BB-level where the 
choice of methods for solving Eq. 2 is autonomous. 

1.5. BLISS Iterative Procedure 

The retrieval of data from SRS (Eq. 3) in a non- 
linear system, is burdened with an error, t* = YV - 
Y A V whose control requires iteration between Eq. 2 
and 5 involving move limits incorporated in QL and 
QU in Eq. 5d. These move limits are being adjusted in 
each iteration and, occasionally, they may push against 
the SRS bounds. That may require the addition of new 
points at which to solve Eq. 2, and refitting of the 
SRS. 

The iterative procedure may also include at its con- 
clusion a retrieval of the optimal XIoc whose elements 
might have been stored in form of SRS. If the volume 
of such storage is prohibitively large, one may choose 
to regenerate the XIoc instead by repeating Eq. 2 w hile 
substituting the latest optimal JXsh|Y*|Y A jw} . How- 
ever, such regeneration is tantamount to a partial execu- 
tion of one additional iteration of the entire BLISS 
procedure, therefore, some of the consistency between 
Eq. 2 and 5 is lost resulting in an additional error. 


i:\amplc: Displacements normal to a wing may be represented by a 
function in the chordwise and spanwise coordinates. Making that 
function linear chordwise and cubic spanwise reduces the number 
of parameters p to only 5. 
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The overall, step-by-step recipe for BLISS 2000 
optimization algorithm is now written as: 

0. Start 

1 . Initialize Xsh and w. and the U and L bounds 

2. Solve equations 

a. h = 0 

b. c = 0 

i.e., perform an initial System Analysis to obtain 
starting values of V* to be considered in 
establishing the RS bounds in the next step. 

3. Approximate model development for each BB. This 
step may be done simultaneously for all the BB's. 

a. Reduce dimensionality of the Q space per Eq. 4. 

b. Disperse by a suitable DOE technique a 
minimum number of points required to define 
an RS approximation in the BB Q-space, 
bounded by QL and QU. 

c. Solve Eq. 2 in subspace U at the above points 
in space Q. This may be done simultaneously 
for all the points within Q. 

d. Fit a Sheaf of Response Surfaces (SRS) to the 
results of #3b. 

e. Verify quality SRS by random sampling. If 
needed, add new points and discard old points, 
and use least square fit (or an equivalent 
technique) to improve the SRS quality. 

f. After each system optimization (#4 below): 
shift, extend or shrink the intervals (QL,QU) to 
avoid excursions beyond the SRS bounds and to 
maintain the approximation quality. 

4. Solve system optimization Eq. 5 in space Q, 
accessing the SRS data per Eq. 3. 

5. Check the termination criteria: Exit, or Repeat from 
#3 using SRS already available or updated per #3e 
and f. 

6. Retrieve the optimal Xloc 

7. Stop 

One execution of the procedure from #3 through 5 
is called the major iteration or the cycle. 

The B11SS Rationale. Theoretical Justification, 
and Computational Effort 

Having the BLISS 2000 algorithm defined, atten- 
tion now turns to its underlying rationale, theoretical 
justification, and the computational cost of the method. 

2,1, Qualitative Discussion of Rationale 

The system objective F is computed as an element 
of Y A in one of the BB's. The Y A from any BB influ- 
ences F because of the couplings c. However, one does 
not initially know the strength or even the sign of these 
influences. The subpotimization task for a BB is, then, 
to develop a wide choice of the BB designs . each 


having a different set of outputs Y A and each Being 
feasible with respect to the local constraints h and g. 

That task is accomplished by using the BB opti- 
mization in Eq. 2 merely as a tool to control the BB 
output Y A , the weighting coefficients w in Eq. 2b 
acting as parameters of that control. Replacement of the 
set of the original vector of the system-level design 
variables ‘Xsh | Y*| with the vector {Xsh | Y* | w{, 
that includes w as additional design variables, results 
in a range of the BB designs, {Xloc | V A {, at each 
point in "the {Xsh | Y*| space, all the BB designs be- 
ing feasible. Availability of a wide choice of these de- 
signs in form of SRS enables the system-level optimi- 
zation in Eq. 5 to find a set of the BB designs that arc 
compatible with respect to c = 0, and generate outputs 
Y A whose collective influence minimizes the objective 
Eq. 5b. 

The above approach to BB suboptimization rests 
on the concept of optimizing a component to attain a 
desired response. That concept may be illustrated by a 
simple example of a cantilever beam of length L and a 
square cross-section of the side length a , loaded by a 
force P at the tip. Consider the beam volume V = a“ L, 
the tip displacement d - P LV(3 E I), where I - a 4 /12, 
and the bending stress at the root s - 6 P L/(a ), as 
three outputs that define the beam response. Let's now 
choose arbitrarily tw'o items out of these three outputs, 
V and d to show hew their control may be achieved 
using a weighting coefficient w inserted in the com- 
posite objective function in the optimization formulated 
as 

GIVEN: constants of E - Young's modulus; 

L - beam length; and s a = allowable 
stress; 

FIND: dimension a 

MINIMIZE: f = V + wd 
SATISFY: s <= s a 

Substituting the formulas for V and d into f, and 
setting dfVdw - 0 yields 

a t - max(Ci w (1 (6PL/s a u )); 

Vo = C: w" do = Cl W' 2 M ; Vo/d» = w M " 

where C, and C; are constants comprising parameters 
P, L, E. and s a . This solution shows how w controls 
the ratio of the two response outputs. 

One should acknowledge that a similar control of 
the outputs may be accomplished by other means. For 
example, the Collaborative Optimization method (e.g.. 
Refs. 9 and 1 1 ) achieves that control by imposing tar- 
gets Y 1 on the output variables Y A and minimizing the 
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discrepancy t— T (Y ! - Y A ) between Y 1 and Y A as the 
objective Junction in Eq. 2 to generate the optimal Y A . 
Unlike BLISS 2000, CO includes the Xsh inputs as 
des ign variables in U in Eq. 2 together with Xloc, and 
imposes the targets on Xsh to generate optimal Xsh in 
addition to optimal Xloc. The vv variables are absent in 
the CO formulation as their role is assigned, in effect, 
to the target variables. In both the CO version reported 
in Ret. 1 1 and in BLISS 2000 herein, the optimal re- 
sults are stored tor the use in the system-level optimiza- 
tion in form of RS. 

The formulation with targets Y 1 in place of w ap- 
pears to be a viable alternative in BLISS. It could be 
combined with stating the objective function in Eq. 2b 
in a Goal Programing style- 2 as f = J (d 4 d ) 
where d and d measure the overachievement or 
underachievement of the targets. Examination of merits 
of the above alternatives is beyond the scope of this 
paper. 

2.2. Theoretical Justification 

The following chain reasoning show's that a solu- 
tion of the set of Eq. 2 and 5 satisfies the original prob- 
lem described by Eq. I: 

1 . Solution of Eq. 2 satisfies constraints h - 0 and 
g <- 0 present in Eq. I . 

2. Solution of Eq. 5 satisfies constraints c — 0 present 
in Eq. I. 

3. Solution of Eq. 2 and the subsequent generation of 
SRS, Eq. 3, renders Y A = Y A “«( {Xsh | Y* | w(), 
and Xloc = Xloc( [Xsh | Y* | w}). If the problem is 
convex, the above relations are unique for each 
vector element (single-valued functions). 

4. Owing to the abov e relations, the space of 

(Xsh j Y* ] w } , searched in Eq. 5, maps uniquely 
on the original space of (Xsh j Xloc | Y* | Y A j, 
searched in Eq. I, because Xsh and Y* define both 
spaces, and Xloc and Y A are turned into functions 
of | Xsh | Y* | w | per #3 above. 

5. The existence of the mapping of the space of Eq. 5 
onto the space of Eq. 1 together with satisfaction of 
the constraints g, h, and c, assures that the solution 
of the sequence Eqs. 2, 3, 5 satisfies also Eq. 1, 
subject to the caveats that follow. 

Caveat 1: If the problem of Eq. 1 has local minima 
(a non-convex problem), then the solutions of Eqs. 2, 

3, 5 may not arriv e at the solution of Eq. 1 even when 
starting from the same initial point because the algo- 
rithmic differences may result in different search path 
being traced. 

Caveat 2: If the SRS were error free and wide 
enough to contain the solution to Eq. 5, a single execu- 
tion of the sequence Eq. 2 Eq. 5 would suffice. 


Otherwise, the sequence has to include the SRS updates 
and be iterated to reduce the error e - Y A „ - Y A V 

2.3. Computational Effort 

in BLISS 2000 the SRS data may be generated 
simultaneously and independently at distributed sites. 
Furthermore, this data generation does not require any 
modifications to the existing codes. Each code involved 
in Eq. 2 may simply be replicated over many proces- 
sors and executed with different inputs Q. It is a classic 
case of the so-called course-grain parallelism in comput- 
ing. Thus, there is no additional up-front reprogram- 
ming cost although there is some cost penalty in com- 
munication and bookkeeping. 

The BLISS computational cost as measured by the 
elapsed time is primarily in repetitive execution of 
Eq. 2 that for a large-scale problem overshadows the 
cost of Eq. 5. For a quadratic polynomial RS in the Q 
space of N dimensions the minimum number of points 
is M s 0(N72). However, one may choose to increase 
the number of points to M 2 = M, M, Ml > 1, to 
improve the RS accuracy, in which case the RS may be 
fitted by means of the method of least squares (or an 
equivalent) with the added benefit of smoothing that 
lessens the probability of entrapment in a local mini- 
mum. Furthermore, the procedure of solving the 
sequence Eq. 2-5 requires a problem-dependent number 
ol cycles, NQ, so the total number of points increases 
to 0(NQ, N72). Consequently, the quadratic RS may 
be impractical beyond N = 12 to 20 in applications 
where the elapsed computing time per point is substan- 
tial and all calculations are sequential. Therefore, the 
use of BLISS 2000 in large-scale applications depends 
on the technology of Massively Concurrent and Dis- 
tributed Processing (MCDP) for compressing M opera- 
tions into the elapsed time of one. Such compression 
not only extends radically the RS limits of practicality 
now but also provides for their further relaxation with 
the progress of the MCDP technology. 

Another way to circumvent the "curse of dimen- 
sionality^ is to settle for a linear form of RS that makes 
M 2 O(N). Then, the iterative process of solving the 
sequence of Eq. 2-5 becomes similar to the sequential 
linear programming (or its more refined version of the 
sequential approximate programming 22 ). Many years of 
experience of the above practice with diverse applica- 
tions indicates that the number of times, NA, the ap- 
proximations must be refreshed in these sequential 
processes depends on the class of application but stays 
independent of N and reasonably small within that 
class. For instance 23 - 4 , in structural optimization 
dominated by membrane stresses, typically. NA = 5 
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to 10, and it rises to 20 to 30 where bending is 
prevalent- 4 . The total number of points to be evaluated 
in such a process is of ()(NA, N). A recent numerical 
evidence attesting effectiveness of the linear RS in a 
sequential approximate optimization is provided in 
Refs. 25 and 26 where the linear RS is used in an 
MCDP environment, and in Ref. 1 1 where it is an 
enhancement in CO. 

Consequently, substituting NA and NQ expected 
for a particular BB and checking whether 0(NA N) < 
Q(NQ, N : /2) is true constitute an important factor in 
the choice between the quadratic or linear RS for that 
particular BB. 

Another factor is an organizational one. The RS 
refreshing operations (Eq. 2) need to be completed at 
the outset of each cycle so a degree of coordination 
among the groups performing these operations is neces- 
sary. The cost of that coordination adds to the cost of 
the entire procedure. On that score, the higher order RS 
are superior because NQ is likely to be smaller than 
NA, hence, the coordination is needed less otten. 

3. A Numerical Example 

A Supersonic Business Jet (SBJ) shown in Fig. 1 
is a numerical example for the BLISS 2000 algorithm. 
It the same case that was used in Refs. 10, 19, and 20 
in the overall layout, but sufficiently different in detail 
to make the results incomparable. Ref. 27 completely 
documents the present case and its relation to the previ- 
ous one. The model of the aircraft system is reduced to 
four BB's shown in Fig. 3 which shows also their data 
exchanges in a standard format of the data dependence 
matrix, also known as the n-square format. That format 
depicts the modules strung along the diagonal. Each 
module accepts input vertically from above or below, 
and outputs horizontally, left or right. 

The BB-structures is a plate representation of the 
wingbox, connected to a rigid beam model of the fuse- 
lage to which a horizontal tail is attached. The wingbox 
strength (stress and local buckling constraints) and 
stiffness accounts for the wing sandwich covers. The 
wingbox volume is the fuel tank volume. The BB- 
aerodynamics calculates lift distribution sensitive to the 
wing trapezoidal geometry' changes and to the local 
variations of the angle of attack generated by an aerody- 
namic twist and by the elastic deformation. This BB 
includes also computation of the drag accounting for 
the wave drag, and accounts tor the trim constraints by 
performing a trim analysis involving the tail volume. 
The BB-power simulates the propulsion by interpola- 
tion of a look-up table that contains data on the thrust 


and specific fuel consumption as functions of the Mach 
number and altitude. Finally, the BB-performance 
computes the flight range by the Breguet equation. 

The system objective is the range under constant 
TOGW. The system design variables govern the geome- 
try and include the Mach number and altitude. The 
local variables are: in structures, the cross-sectional 
dimensions of the sandwich wing covers at several loca- 
tions over the wing; in aerodynamics, the horizontal 
tail sweep angle, wing and horizontal tail location co- 
ordinates; and in Power the throttle setting. Table 1 
provides more information for the input, output, and 
the RS statistics for the BB s. The field variables of 
aerodynamic loads and the wing twist were condensed 
to their distribution parameters. 

The n-square diagram in Fig. 4 corresponds to the 
one in Fig. 3, except that it shows the system level 
optimizer as the sole input source tor all BB s. Each 
BB in that diagram is now a SRS representation of that 
BB suboptimization. The c-constraints occur at the data 
channel intersections marked by circles. 

The design points for construction of RS are placed 
using the D-optimal technique from the Design of 
Experiments methodology. The RS are quadratic poly- 
nomials. As the procedure progresses, they are being 
periodically updated in order to keep them centered in 
the Q space around the latest optimal solution, to pre- 
vent the optimal solution to remain lodged against the 
RS boundary, and to reduce the RS span as the proce- 
dure homes on the system optimum. 

The BB analysis codes are in FORTRAN, while 
MATLAB Optimization Toolbox is the tool for opti- 
mization at the BB level and the system level. The 
main program that executes the entire procedure is also 
written in MATLAB. 

In this example, BLISS 2000 manipulates the sys- 
tem and local design variables so as to adjust the wing 
structural weight, engine weight, aerodynamic drag, 
fuel volume, and fuel consumption so as to maximize 
the range. 

The results are very voluminous- 7 so only a small 
but representative sample is displayed herein in Table 2 
that shows how the objective and a few' local and sys- 
tem variables were changing over the BLISS cycles. In 
this case, the result variations diminish to the level 
commensurate with the accuracy of the analysis typi- 
cally after 6 or 8 cycles, hence the number of cycles 
is set to ten. The data are shown normalized by their 
values obtained in Iteration 10; the latter are also pre- 
sented in the right-most column. The BB internal op- 
timization histograms are not shown because the local 


7 

American Institute of Aeronautics and Astronautics 



optimization performance has no hearing on the overall 
convergence of the method owing to the BB autonomy. 

fable 2 includes the objective function of the flight 
range. At a cycle number, the table shows the range 
predicted by RS, Rrs, used in that cycle and the range 
actually computed. Ran, at the outset of the next cycle. 
The two values converge into a complete agreement in 
five cycles. The convergence is oscillatory partially due 
to the initial differences of the RS-predicted and actual 
values and partially due to the search of a compromise 
among the multitude of trade-offs at two levels account- 
ing for constraints. 

The method transforms the initially infeasible de- 
sign to one where all the constraints are satisfied and 
the objective is improved by 47 %. That improvement 
value is significant merely as a metric of performance of 
the method in this particular test but the range final 
value must not be interpreted as realistic. The analysis 
scope and fidelity, and completeness of the set of vari- 
ables and constraints would have to be much higher for 
such interpretation. 

4. Concluding Remarks 

BLISS 2000 method is intended specifically for 
support of design projects that require collaboration of 
autonomous, distributed specialty groups and exploita- 
tion of the Massively Concurrent and Distributed Proc- 
essing (MCDP) technology. 

The method decomposes the overall system opti- 
mization into subotimizations associated with the spe- 
cialty groups. Response Surfaces approximate the 
suboptimization optimal results in the coordinating, 
system-level optimization that, provably. leads to the 
same solution, subject to the approximation errors, that 
would have been obtained without decomposition. An 
example of a simplified, conceptual design of an aircraft 
illustrates the method. 

The proof-of-the-concept testing reported herein 
showed BLISS performing as intended. Preliminary 
assessment of the method merits and demerits follows. 

The BLISS 2000 performance depend s on the qual- 
ity of approximation of the subsystem optimization 
results as functions of the system-level design variables 
and coupling variables. The number of times these ap- 
proximation need to be updated and the number of new 7 
design points in these updates are problem dependent. 
Assessment of the approximation errors is a factor in 
choosing between the linear or higher order approxima- 
tions for the individual modules, and in deciding on 
the frequency of the approximation updates. 


BLISS 2000 tends to generate a large amount of 
computing and voluminous intermediate data. On the 
other hand it simplifies the entire procedure conceptu- 
ally to reduce the human effort and time. This appears 
to be a cost-effective trade-off, considering that comput- 
ing gets less expensive all the time while the labor 
costs grow. The elapsed time consumed by the BL1SS- 
generated computation is compressed by the MCDP 
technology. That technology has established itself the 
future of computing in science and engineering and 
BLISS being intrinsically coarse-grain parallel is posed 
to exploit that technology and grow 7 with it. 

The BLISS method performs best in applications 
that decompose into subsystems w'ith numerous local 
design variables and constraints leaving a relatively 
small number of design variables at the system level. A 
nearly complete autonomy of the operations in the sub- 
systems is the key BLISS advantage that enables full 
utilization of the local know ledge and control over the 
budget and task time. 

The system level-optimization benefits from nearly 
instantaneous response of the approximations, therefore, 
it does not call for any particularly high sophistication 
or efficiency of the search technique. 

Further research is recommended to determine how 
much BLISS could reduce the elapsed project time, and 
whether it would benefit from the alternative formula- 
tions of the module-level optimization, use of the linear 
instead of quadratic form of the Response Surfaces, and 
approximations other than Response Surfaces. It is 
through generation of experience in large-scale practical 
applications that guidelines for further gains in effi- 
ciency and accuracy will develop. 
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Table 1. Variables Z, Y*, Y A and dimensionality data for RS 




Input Variables 


Number 
of Points 

BB 

Output Y a 

Z 

Y* 

W 

Number 
of Inputs 

1. Structures 

1. W T 

2. W ( . 

3. © 

t/c, AR W , 
Sreb 

ARht, ^ 

l,w e 

I 

R 

11 

78 

2. Aero 

4. L 

5. D 

6. L/D 

t/c, h, M, AR W , 
S RI ;p, Sht, 
ARht, ^ 

w T ,0 

ESF 

w2(l) 

w2(3) 



3. Power 

1. SCF 

2. W E 

3. ESF 

h, M 

D 

w3( 1 ) 
w3(3) 

5 

1 

4. Performance 

I . Range 

h, M 

w T , 

W F , 

L/D, 

SFC 


6 



Table 2. A sample of data changing over 10 BLISS cycles 


Cycle 

1 0 1 1 

□T 

nr 

4 

5 

1 — 7 

1 7 

r^r 

FT 

r 10 




Normalized b\ 

Value from Cycle 10 





M 

1 .08 

1.01 

1.07 

1.10 

1.05 

1.03 

1.00 

0.99 

1.00 

0.99 

1.00 

TOG W (lb) 

1.06 

LI 1 

1.19 

1.12 

1.05 

1.02 

1.02 

1.01 

1.01 

1.00 

1.00 

Sref ( ft A 2) 

0.69 

0.35 

0.77 

0.83 

0.77 

0.72 

0.85 

0.86 

0.93 

0.98 

1.00 

ESF 

1.94 

1.22 

0.90 

1.05 

1.10 

1.03 

0.97 

1.00 

1.03 

1.01 

1.00 

t/c 

2.16 

2.89 

2.89 

1.60 

1.69 

1.55 

1.31 

1.26 

1.12 

1.02 

1.00 

tsl(in) d 

1.04 

1.03 

1.05 

1.07 

1.04 

1.00 

0.98 

1.00 

0.99 

0.99 

1.00 

T (%) 

1.70 

1 .03 

1.52 

1.22 

1.10 

1.07 

1.01 

0.98 

1.00 

1.00 

1.00 

Lht (%MAC) 

0.71 

0.61 

1.00 

1.00 

0.60 

0.76 

1.00 

1.00 

1.00 

1.00 

LOO 

w 1 

0.00 

0.00 

0.48 

1.19 

1.90 

1.67 

2.14 

1.43 

1.10 

1.00 

1.00 

vv3 

0.00 

1.79 

1.61 

1.02 

1.16 

0.89 

1.25 

1.02 

1.00 

1.02 

LOO 

Rrs 


0.77 

0.93 

1.20 

0.76 

0.91 

0.97 

0.99 

0.98 

0.99 

1 .00 

Ran(nm) 

a i • a >• 


0.68 

0.91 

1.07 

0.77 

0.91 

0.97 

0.99 

0.98 

0.99 

LOO 


Value from | 
Cycle 10 


1 .663 
31342 
578 
0.517 
0.035 
3.840 
20.6 
350 
0.210 
0.560 
5247 
5247 


sandwich caliper thickness of the win^ upper cover inboard segment 
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Fig. 1. A generic Supersonic Business Jet as a test case 
in BLISS. 



Fig. 2. An example of a few principal disciplines in an 
aircraft treated as a system. 



Fig. 3. Test case coupling data How. 



Fig. 4. System optimizer is a source ot all BB inputs Z 
and Y*. BB’s return Y A to the optimizer (Table 1 
shows Y \ Z, and Y*). 
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